home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / MATH / FSULTRA1.ZIP / ULTRX_TP.ASM < prev    next >
Assembly Source File  |  1992-06-18  |  3KB  |  136 lines
  1. comment ! 
  2. FSU - ULTRA    The greatest random number generator that ever was
  3.         or ever will be.  Way beyond Super-Duper.
  4.         (Just kidding, but we think its a good one.)
  5.  
  6. Authors:    Arif Zaman (arif@stat.fsu.edu) and
  7.         George Marsaglia (geo@stat.fsu.edu).
  8.  
  9. Date:        27 May 1992
  10.  
  11. Version:    1.05
  12.  
  13. Copyright:    To obtain permission to incorporate this program into
  14.         any commercial product, please contact the authors at
  15.         the e-mail address given above or at
  16.  
  17.         Department of Statistics and
  18.         Supercomputer Computations Research Institute
  19.         Florida State University
  20.         Tallahassee, FL 32306.
  21.  
  22. See Also:    README        for a brief description
  23.         ULTRA.DOC    for a detailed description
  24.  
  25. -----------------------------------------------------------------------
  26. ! 
  27. ;
  28. ; File: ULTRX_TP.ASM  (80386 Turbo Pascal calling conventions)
  29. ;
  30.     DOSSEG
  31.    .MODEL TPASCAL
  32.  
  33. ;===== B: GLOBAL FUNCTIONS, SUBROUTINES and VARIABLES =================
  34. ;
  35.     public  i31bit, i15bit, i7bit,  i1bit,  uni,    duni
  36.     public  i32bit, i16bit, i8bit,  rinit,  vni,    dvni
  37.  
  38. ;===== D. MACRO DEFINITIONS ===========================================
  39. ;
  40. ; RinitProcStart should take two 32-bit arguments conx and shrx
  41. ;   and place them in eax and ebx.
  42. ;   es and ds should both point to the data segment.
  43. ;   conx must be odd (so we or with 1 just to make sure).
  44. ; FillProc
  45. ;   DS is already the data segment. ES should also be made the same.
  46.  
  47. EnterProcedure   macro
  48. endm
  49.  
  50. ExitProcedure    macro
  51.     ret
  52. endm
  53.  
  54. Enter2arg   macro
  55.   arg conx:dword, shrx:dword
  56.   EnterProcedure
  57.     mov ax,ds
  58.     mov es,ax
  59.     mov eax,conx
  60.     shl eax,1
  61.     or  al,1
  62.     mov ebx,shrx
  63. endm
  64.  
  65. Exit2arg     macro
  66.   ExitProcedure
  67. endm
  68.  
  69. EnterFill    macro
  70.     mov    ax,ds
  71.     mov    es,ax
  72. endm
  73.  
  74. ExitFill    macro
  75.     ret
  76. endm
  77.  
  78. DwordFn macro
  79. endm
  80.  
  81. WordFn  macro
  82. endm
  83.  
  84. ByteFn  macro
  85. endm
  86.  
  87. RealFn  macro
  88. endm
  89.  
  90. DoubleFn macro
  91. endm
  92.  
  93. ;===== DATA ========================================================
  94.  
  95.   N equ 37              ; The number of 32 bit words in the table
  96.   N2 equ 24             ;  x(i) = x(i-N) - x(i-N2)
  97.  
  98. counterpoint    STRUC
  99.   c  dw 0               ; counter
  100.   p  dw ?               ; offset for pointer
  101.      dw ?               ; segment for pointer
  102.   x  dd N dup (?)    ; data
  103. ENDS
  104.  
  105.   extrn swbstate            ; In TP, the shared DATA is
  106.   cps = size counterpoint
  107.   cp1 = cps - 4*N + 32
  108.   swb32 = counterpoint ptr swbstate    ; in the TPU, so we define a
  109.   swb16 = counterpoint ptr swb32+cps    ; large array there, and refer
  110.   swb8  = counterpoint ptr swb16+cps    ; to it here. This ensures that
  111.   swb1  = counterpoint ptr swb8 +cps    ; the swbstate is contiguous.
  112.   swbseed = dword ptr swb1     + cp1
  113.   flags   = byte  ptr swbseed  + 4*N
  114.   congx   = dword ptr flags    + 1
  115.  
  116.   extrn tmpq  : qword
  117.   extrn neg31 : word
  118.   extrn neg63 : word
  119.  
  120.   conglo  = word ptr congx
  121.   conghi  = word ptr congx+2
  122.   tmpw1   = word ptr tmpq
  123.   tmpw2   = word ptr tmpq+2
  124.   tmpw3   = word ptr tmpq+4
  125.   tmpw4   = word ptr tmpq+6
  126.   tmpdlo  = dword ptr tmpq
  127.   tmpdhi  = dword ptr tmpq+4
  128.   shrglo  = tmpw1
  129.   shrghi  = tmpw2
  130.   shrgx   = tmpdlo
  131.  
  132. .CODE
  133.   INCLUDE ULTRXCOD.INC
  134. END
  135.  
  136.